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Generalized  Korteweg  -  de  Vries  equation  (GKdV).  This  equation  is 
written  as  ^t+^  ^x+'^xxx  -0  where  p  is  a  positive  integer. 

(p  «  1  gives  the  KdV  equation).  -^The  convergence  estimates  obtained 
for  the  KdV  equation  have  been  generalized,  under  the  assumption  that 
the  solution  u  is  sufficiently  regular.  For  p  4,  it  is  not  known 
whether  a  global  smooth  solution  exists  corresponding  to  smooth 
initial  data.  It  is  in  fact  conjectured  that  for  these  cases,  the 
solution  may  develop  a  singularity  in  finite  time.  A  code  that  uses  a 
spatially  and  temporally  adaptive  strategy  has  been  implemented.  We 
are  currently  investigating  the  stability  of  solitary  type  solutions. 
As  conjectured,  these  solutions  are  highly  unstable  for  initial 
amplitudes  larger  than  one.-^  With  the  code,  we  are  able  to  track  the 
solution  to  the  point  where  it  develops  a  pulse-like  shape  with  an 
amplitude  of  2.25  x  10^  supported  over  an  interval  of  width  10-15 
(See  attached  graphs).  These  experiments  will  be  reported  in  (3J. 
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1.  The  Stokes  prohleiii.  Nonconforming  finite  element  approximations  to  solutions  of 
the  Stokes  equations  are  constnicted  in  {IJ.  Optimtd  rates  of  convergence  are  proved 
for  the  velocity  and  pressure  approximations,  in  the  usual  Sobolev  norms  of  positive 
as  well  as  negative  index.  For  the  pressure,  C"  piecewise  polynomial  functions  are 
used.  Tlte  class  of  vector  fields  used  to  approxunate  the  velocity  field  have  piecewise 
polynomial  components,  discontinuous  across  intereleinent  Imundaries.  On  each 
triangle,  these  vector  fields  satisfy  the  incompressibility  condition  pointwise.  It  is  shown 
that  these  piecewise  solenoidal  vector  fields  possess  optimal  approximation  properties 
to  smooth  solenoidal  vector  fields  on  domains  with  curved  boundaries.  Bases  for  these 
approximating  spaces  can  be  constnif'ted  very  easily,  and  bases  for  higher  order  spaces 
can  be  constructed  by  adding  suitable  tenns  to  those  for  spaces  of  lower  order. 


2.  Kortewog  -  do  ^Vic.s  pf|natioii.  (KdV).  In  [2j,  we  show  that  the  well-known  order 
reduction  phenomenon  affectmg  Implicit  Runge-Kutta  methods  does  not  occur  when 
approximating  periodic  solutions  of  the  KdV  equation,  i.e.  that  the  rates  of  convergence 
are  those  that  one  expects  in  the  case  of  nonstiff  systems  of  ODE’s.  For  example,  if  a 
q-stage  IRK  method  of  Gauss-Legendre  type  is  used,  then  the  rate  of  convergence  is 
2q.  Interestingly,  our  estimates  may  be  viewed  as  superconvergent  if  the  IRK  method 
is  viewed  as  a  special  case  of  a  continuous  finite  element  in  time  fonnulation. 


3.  Geiiet  nllzed  Kortewog  -  cle  Vries  equation  (GKdV).  This  equation  is  written  as 

II,  I  ll''ll  I-  II,  rr  —  (I 

where  y»  is  a  positive  integer.  (y»  =  1  gives  the  KdV  equation).  The  convergence 
estimates  obtained  for  the  KdV  ec|uation  have  been  generalized,  under  the  assumption 
that  the  solution  n  is  sufficiently  regular.  For  />  -  1,  it  is  not  known  whether  a  global 
smooth  solution  exists  corresponding  to  smooth  initial  data.  It  is  in  fact  conjectured 
that  for  these  ca.ses,  the  solution  may  develop  a  singularity  in  finite  tune.  A  code 
that  uses  a  spatially  and  temporally  adaptive  strategy  has  been  implemented.  We  are 
I  currently  investigating  the  stability  of  solitary  type  solutions.  As  conjectured,  these 
solutions  are  highly  unstable  for  initial  amplitudes  larger  than  one.  With  the  code,  we 
me  able  to  track  the  solution  to  the  point  where  it  develops  a  pulse-like  shape  with 
an  amplitude  of  2.'2~)  x  10’  supported  over  an  interval  of  width  10  '  (See  attached 
graphs).  These  experiments  will  be  reported  in  [3]. 


4.  Noiiltnear  (C'ubir)  Soiirodinger  f»nia(ioii  (NLS).  We  consider  the  initial¬ 
boundary  valur  problem 
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here  U  is  a  bounded  domaitt  in  R"  and  \  is  a  real  parameter.  In  [4],  Tire  solution  n  is 
approximated  by  tlie  Galerkin  method  in  space  and  two  second-order  accurate  temporal 
discretizations,  fn  15],  high  ord^r  temporal  discretizations  generated  by  IRK  methods 
are  analyzed.  For  general  domains,  a  rate  of  convergence  of  </  f  ■]  is  shown  to  hold, 
wliile  for  the  case  of  rectangular  domains  with  sides  perpendicular  to  the  cr>ordinate 
axes,  (including  the  important  one  dimensional  case)  the  full  rate  of  '2<i  comes  into 
force.  The  computational  work  consists  in  solving,  at  each  time  level,  a  set  of  </  linear 
systems.  These  matrices,  albeit  different  from  one  another,  are  independent  from 
time,  re.sulting  in  a  very  efficient  scheme.  In  addition,  these  </  systems  may  be  .solved 
independently  on  a  “pardlel”  computer  with  at  least  </  processors. 

5.  P.'u-olh'l  mot  Inxl.s.  The  inherent  parallel  features  of  IRK  methods  are  exploited, 
resulting  in  schemes  that  give  high  accuracy  while  requiring  execution  times  similar 
to  those  of  the  lowest  order  methods.  Such  schemes  for  the  KdV  and  NLS  equations 
have  already  been  implemented  and  tested  on  the  CRAY  and  the  Hypercube,  two 
machines  that  typify  popular  but  different  philosopliies  in  parallel  computing.  Speedup 
measurements  are  under  way,  as  well  as  the  task  of  gathering  data  regarding  the  relative 
efficiency  of  high-order  nteilunls  versus  those  of  low  order.  Tlie  conclusions  of  the.se 
experiments  will  be  reported  in  the  near  future. 
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